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Abstract. A local algorithm is a distributed algorithm where 
each node must operate solely based on the information that was 
available at system startup within a constant-size neighbourhood 
of the node. We study the applicabiUty of local algorithms to 
max-min LPs where the objective is to maximise min^ CkvXy 
subject to aiyXy < 1 for each i and x^^ > for each v. Here 
Ckv > 0, ttii, > 0, and the support sets Vi ~ {v : > 0}, 
Vfc — {v : Ckv > 0} , Iv — {i ■ CLiv > 0} and = {k : Ckv > 0} have 
bounded size. In the distributed setting, each agent v is responsible 
for choosing the value of x^ , and the communication network is a 
hypergraph 'K where the sets Vk and Vi constitute the hyperedges. 
We present inapproximability results for a wide range of struc- 
tural assumptions; for example, even if and \Vk\ are bounded 
by some constants larger than 2, there is no local approximation 
scheme. To contrast the negative results, we present a local ap- 
proximation algorithm which achieves good approximation ratios 
if we can bound the relative growth of the vertex neighbourhoods 
in J{. 



1. Introduction 

We study the limits of what can and what cannot be achieved by 
local algorithms |13|. We focus on the approximability of a certain 
class of linear optimisation problems, which generalises beyond widely 
studied packing LPs; the emphasis is on deterministic algorithms and 
worst-case analysis. 

1.1. Local algorithms. A local algorithm is a distributed algorithm 
where each node must operate solely based on the information that 
was available at system startup within a constant-size neighbourhood 
of the node. We focus on problems where the size of the input per 
node is bounded by a constant; in such problems, local algorithms 
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provide an extreme form of scalability: the communication, space and 
time complexity of a local algorithm is constant per node, and a local 
algorithm scales to an arbitrarily large or even infinite network. 

The study of local algorithms has several uses beyond providing 
highly scalable distributed algorithms. The existence of a local al- 
gorithm shows that the function can be computed by bounded-fan-in, 
constant-depth Boolean circuits; we can say that the function is in the 
class NC°. A local algorithm is also an efficient centralised algorithm: 
the time complexity of the centralised algorithm is linear in the num- 
ber of nodes; furthermore, due to spatial locality in memory accesses, 
we may be able to achieve a low I/O complexity in the external mem- 



ory 18 model of computation. In certain problems, a local approxi- 
mation algorithm can be used to construct a suhlinear time algorithm 
which approximates the size of the optimal solution, assuming that we 
tolerate an additive error and some probability of failure jlG]. A local 
algorithm can be turned into an efficient self- stabilising algorithm jsj; 
the time to stabilise is constant jl] . Finally, the existence and nonexis- 
tence of local algorithms gives us insight into the algorithmic value of 
information in distributed decision-making [14j. 

1.2. Max-min packing problem. In this section, we define the opti- 
misation problem that we study in this work. Let V , I and K be index 
sets with J n -ft" = 0; we say that each u e V" is an agent, each k & K is a. 
beneficiary party, and each z e / is a resource {constraint). We assume 
that one unit of activity by v benefits the party k by Ckv > units and 
consumes Oj^, > units of the resource i] the objective is to set the 
activities to provide a fair share of benefit for each party. In notation, 
assuming that the activity of agent v is x„ units, the objective is to 

maximise u = min } CkvXy 

(1) subject to Q-iv^v < 1 for each i G /, 

Xt, > for each v ^V. 

Throughout this work we assume that the support sets defined for 
alH e /, /c e K, and v ^ V hy Vi = {v ^ V : ai^ > 0}, 14 = {t; G 
V : > 0}, J„ = {i G / : a™ > 0}, and = {k ^ K : Ckv > 0} 
have bounded size. That is, we consider only instances of Q such that 
\Iv\ < ^v, \Kv\ < Ay, \Vi\ < Aj and \Vk\ < A^ for some constants 
Ay, Ay, Aj and A]^. To avoid uninteresting degenerate cases, we 
furthermore assume that J„, Vi and 14 are nonempty. 
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1.3. LP formulation. If the sets V, I and K are finite, the problem 
can be represented using matrix notation. Let A be the nonnegative 
|/| X \V\ matrix where the entry at row i, column v is aj^,; define C 
analogously. We write for the row i of A and Ck for the row k of 
C. Let column vector of length \V\. The goal is to maximise 
LJ = miufcgj^ CkX subject to Ax < 1 and x > 0. 

In the special case \K\ = 1, this is the widely studied fractional 
packing problem: maximise cx subject to Ax < 1 and x > 0. This 
simple linear program (LP) has nonnegative coefficients in c and A. We 
refer to a problem of this form as a packing LP; the dual is a covering 
LP. Naturally the case of any finite K can also be written as a linear 
program, but the constraint matrix is no longer nonnegative: maximise 
u subject to Ax < 1, ul — Cx < and x > 0. 

1.4. Distributed setting. We construct the hypergraph CK = {V,8,) 
where the hyperedges are £ = {Vi : i E 1} U {Vk : k G K}. This is the 
communication graph in our distributed optimisation problem. The 
variable x^ is controlled by the agent v E V, and two agents u,v E V 
can communicate directly with each other if they are adjacent in "K. 
We write dji{u, v) for the shortest-path distance between u and v in 
%. The agents are cooperating, not selfish; the difficulty arises from 
the fact that the agents have to make decisions based on incomplete 
information. 

Initially, each agent v eV knows only the following local informa- 
tion: the identity of its neighbours in the graph J{; the sets 1^ and K^] 
the values for each i E ly] and the values Ckv for each k G K^. That 
is, V knows with whom it is competing on which resources, and with 
whom it is working together to benefit which parties. 

When we compare the present work with previous work, we often 
mention the special case \K\ = 1, as this corresponds to the widely 
studied packing LP. However, in this case the size of Vk is not bounded 
by a constant A]^: we have Vk = V for the sole k E K. Therefore 
we introduce a restricted variant of the distributed setting, which we 
call collaboration-oblivious. In this variant, the hyperedges are £ = 
{Vi : i G /}. Whenever we study related work on the packing LP, we 
focus on the collaboration-oblivious setting. 

1.5. Local setting. We are interested in solving the problem ([T| by 
using a local algorithm. Let r = 1,2,... be the local horizon of the 
algorithm; this is a constant which does not depend on the particular 
problem instance at hand. Let B^jiiy^r) = {u eV : d'xiu^v) < r} be 
the set of nodes which have distance at most r to the node v in IK. 
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The agent v must choose the value based on the information that is 
initially available in the agents B'ji{v,r). 

We focus on the case where the size of the input is constant per node. 
The elements and Ckv are represented at some finite precision. Fur- 
thermore, we assume that the nodes have constant-size locally unique 
identifiers; i.e., any node can be identified uniquely within the local 
horizon. 

1.6. Approximation. A local algorithm has the approximation ratio 
a for some a > 1 if the decisions x^ are a feasible solution and the 
value uj is within a factor a of the global optimum. A family of local 
algorithms is a local approximation scheme if we can achieve any a > 1 
by choosing a large enough local horizon r. 

1.7. Contributions. In Section |4] we show that while a simple algo- 
rithm achieves the approximation ratio A]^ for ([T]), no local algorithm 
can achieve an approximation ratio less than Ay/2 + 1/2 — 1/(2A^ — 2) 
in the general case. In Section |5] we present a local approximation al- 
gorithm which can achieve an improved approximation ratio if we can 
bound the relative growth of the vertex neighbourhoods in 

2. Applications 

Consider a two-tier sensor network: battery-powered sensor devices 
generate some data; the data is transmitted to a battery-powered relay 
node, which forwards the data to a sink node. The sensor network is 
used to monitor the physical areas K. Let 5* be the set of sensors, and 
let T be the set of relays; choose I = S UT. 

For each sensors device s G 5*, there may be multiple relays t E T 
which are within the reach of the radio of s; we say that there is a 
wireless link (s, t) from s to t. The set V consists of all such wireless 
links, and the variable X{^s,t) indicates how much data is transmitted 
from s via t to the sink. Transmitting one unit of data on the link 
V = {s, t) E V and forwarding it to the sink consumes the fraction a^^ 
of the energy resources of the sensor s and also the fraction of the 
energy resources of the relay t. 

Let Ckv = 1 for each link v = {s,t) if the sensor s is able to monitor 
the physical area k E K. Now (jlj) captures the following optimisation 
problem: choose the data flows in the sensor network so that we max- 
imise the minimum amount of data that is received from any physical 
area. Equivalently, we can interpret the objective as follows: choose 
data flows such that the lifetime of the network (time until the first 
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sensor or relay runs out of the battery) is maximised, assuming that 
we receive data at the same average rate from each physical area. 

Similar constructions have applications beyond the field of sensor 
networks: consider, for example, the case where each k & K is a major 
customer of an Internet service provider (ISP), each s G S" is a bounded- 
capacity last-mile link between the customer and the ISP, and each 
t G T is a bounded-capacity access router in the ISP's network. 



3. Related work 



Papadimitriou and Yannakakis 15 present the safe algorithm for 
the packing LP. The agent v chooses 

1 

(2) Xv = mm — — . 

«e/„ aiv\Vi\ 

This is a local A]^-approximation algorithm with horizon r = 1. 

Kuhn et al. [o] present a distributed approximation scheme for the 
packing LP and covering LP. The algorithm provides a local approx- 
imation scheme for some families of packing and covering LPs. For 
example, let G {0, 1} for all i, v. Then for each , Ay and a > 1, 
there is a local algorithm with some constant horizon r which achieves 
an a-approximation. Our work shows that such local approximation 
schemes do not exist for ([T]). 

Another distributed approximation scheme by Kuhn et al. [9] forms 
several decompositions of "K into subgraphs, solves the optimisation 
problem optimally for each subgraph, and combines the solutions. How- 
ever, the algorithm is not a local approximation algorithm in the strict 
sense that we use here: to obtain any constant approximation ratio, 
the local horizon must extend (logarithmically) as the number of vari- 
ables increases. Also Bartal et al. j2] present a distributed but not local 
approximation scheme for the packing LP. 

Kuhn and Wattenhofer ^lOj present a family of local, constant-factor 
approximation algorithms of the covering LP that is obtained as an LP 
relaxation of the minimum dominating set problem. Kuhn et al. [T] 
present a local, constant-factor approximation of the packing and cov- 
ering LPs in unit-disk graphs. 

There are few examples of local algorithms which approximate linear 
problems beyond packing and covering LPs. Kuhn et al. {sj study an 
LP relaxation of the fc-fold dominating set problem and obtain a local 
constant-factor approximation for bounded-degree graphs. 

For combinatorial problems, there are both negative [6 ,11 and pos- 
itive |4|[8|[To|[T3|[T7] results on the applicabihty of local algorithms. 
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4. Inapproximability 

Even though the safe algorithm [TEj was presented for the special 
case of li^'l = 1, c = 1, and finite / and V, we note that the safe 
solution X defined by ^ and the optimal solution x* also satisfy 

min Cfc„x* < min ^ Ckv^Jx^ = Aj min ^ CkvXy. 

k&K ^ k&K ^ k(^K ^ 

Therefore we obtain a local approximation algorithm with the approx- 
imation ratio Ay for ([T|. 

One could hope that widening the local horizon beyond r = 1 would 
significantly improve the quality of approximation. In general, this 
is not the case: no matter what constant local horizon r we use, we 
cannot improve the approximation ratio beyond Ay/2. In this section, 
we prove the following theorem. 

Theorem 1. Lei A^ > 2 and A\ > 2 he given. There is no local 
approximation algorithm for ([T]) with the approximation ratio less than 
Ay/2 + 1/2 — 1/(2 Ay — 2). This holds even if we make the following 
restrictions: G {0, 1}, Ay = 1 and Ay = 1. 

We emphasise that the local algorithm could even choose any local 
horizon r depending on the bounds Ay, Ay, Ay and Ay. Never- 
theless, an arbitrarily low approximation ratio cannot be achieved if 
Ay > 3 or Ay > 3. In the case Ay = Ay = 2 the existence of a local 
approximation scheme remains an open question. 

Analogous proof techniques, using constructions based on regular bi- 
partite high-girth graphs, have been applied in previous work to prove 
the local inapproximability of packing and covering LPs |9| and com- 
binatorial problems \6]. 

4.1. Proof outline. Choose any local approximation algorithm /l for 
the problem ([T]). Let r > 1 be the local horizon of A and let a be the 
approximation ratio of A. We derive a lower bound for a by construct- 
ing two instances of ([T]), S and such that certain sets of nodes in 
the two instances have identical radius-r neighbourhoods in both in- 
stances. Consequently, the deterministic local algorithm A must make 
the same choices for these nodes in both instances. The nodes with 
identical views are selected based on the solution of § computed by A, 
which enables us to obtain a lower bound on a by showing that this 
solution is necessarily suboptimal clS 8b solution of S'. 

4.2. Construction of §. We now proceed with the detailed construc- 
tion of the instance §. The constructions used in the proof are illus- 
trated in Figure [1} 
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(a) 




(b) 



Hypcrcdgcs of type II 



A leaf node, 
level 5 



> Hypcrcdges of type I 
The root node, level 




Figure 1. The construction of §, in the case c? = 2, 
D = 3, r = 2, R = 3. (a) A small part of the bipartite, 
72-regular, high-girth graph Q. (b) A complete (2, 3)-ary 
hypertree of height 5, with 72 leaves, (c) The underlying 
hypergraph of S. Grey highlighting indicates the under- 
lying hypergraph of black circles are the variables of 
S' which we set to 1 in Section 14.51 
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Let Aj > 2 and A]^ > 2; without loss of generality we can assume 
that at least one of the inequalities is strict because setting A]^ = 
A^ = 2 in the theorem statement yields the trivial bound a > 1. Let 
d = Ay - 1 and D = A^ - 1. Observe that dD > 1. Let R > r; the 
precise value of R is chosen later and will depend on d, D and a only. 

Let Q be a d^D^~^-Tegu\ai bipartite graph with no cycles consisting 
of less than 4r + 2 edges. (A random regular bipartite graph with suf- 
ficiently many nodes has this property with positive probability [l2].) 
The graph Q provides the template for constructing the hypergraph 
underlying the instance S. 

Before describing the construction, we first introduce some terminol- 
ogy. A complete {d,D)-ary hypertree of height h is defined inductively 
as follows. For h = 0, the hypertree consists of exactly one node and 
no edges; the level of the node is 0. For h > 0, start with a complete 
{d, D)-ary hypertree of height h — 1. For each node v at level h—1, in- 
troduce a new hyperedge and new nodes as follows. If — 1 is even, the 
new hyperedge consists of the node v and d new nodes. U h — 1 is odd, 
the new hyperedge consists of the node v and D new nodes. For future 
reference, call these hyperedges of types I and II, respectively. The 
new nodes have level h in the constructed hypertree. The constructed 
hypertree is a complete {d, D)-ary hypertree of height h. The root of 
the hypertree is the node at level 0, the leaves are the nodes at level 
h. Each level i has either {dDy/"^ or {dD)^^~^'>^'^d nodes depending on 
whether ^ is even or odd, respectively. See Figure [T] for an illustration. 

We now construct the hypergraph underlying S. Denote by Q the 
vertex set of Q. Form a hypergraph "K by taking \Q\ node-disjoint 
copies of a complete [d, D)-aiy hypertree of height 2R — 1. For q E Q, 
denote the copy corresponding to g by Tg. Denote the node set of 7q 
by Tq. For £ = 0, 1, ... , 2R — 1, denote the set of nodes at level i in 7q 
by Tq{i). Denote the set of leaf nodes in 7q by Lq = Tq{2R — 1). 

Observe that the number of leaf nodes in each 7q is equal to the 
degree of every vertex in Q. For each vertex q & Q and each leaf node 
V & Lq, associate with v a unique edge of Q incident with the vertex 
q. Each edge of Q is now associated with exactly two leaf nodes; by 
construction, these leaf nodes always occur in different hypertrees 7q. 
For a leaf v G UgLg, let f{v) be the other leaf associated with the 
same edge of Q. Observe that f{f{v)) = v holds for all v G UqLq] in 
particular, / is a permutation of UqLg. To complete the construction 
of "K, add the hyperedge {v, f{v)} to J{ for each v G UqLq. Call these 
hyperedges type III hyperedges. 

Let us now define the instance S of ([T]) based on the hypergraph !K. 
Let the set of agents V be the node set of J-C. For each hyperedge e of 
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type I, there is a resource i G /; let ai^, = 1 if f G e, otherwise a^^ = 0. 
For each hyperedge e of type II, there is a beneficiary party k & K; let 
Cfe^ = 1/D if f G e, otherwise Ckv = 0. For each hyperedge e of type 
III, there is a beneficiary party k & K; let c^^ = 1 if f G e, otherwise 
Ckv = 0. The locally unique identifiers of the agents can be chosen 
in an arbitrary manner. (This proof applies also if the identifiers are 
globally unique; for example, we can equally well consider the standard 
definition where the identifiers are a permutation of 1, 2, . . . , |V^|.) This 
completes the construction of S. Observe that S has "K as its underlying 
hypergraph. 

4.3. Construction of Next we construct another instance of ([T]), 
called §', by restricting to a part of §. To select the part, we apply 
the algorithm A to the instance §. We do not care what is the optimal 
solution of S; all that matters at this point is the fact that each agent 
V & V must choose some value > 0. In particular, we pay attention 
to the values x„ at the leaf nodes v G U^Lg. 
For all g G Q, let 



For all P Q, let S{P) = ^ggp5(g). Because / is a permutation of 
UqLq with /(/(f)) = V for all v G UqLg, we have S{Q) = 0. Thus, there 
exists a p G Q with S{p) > 0. 

The instance S' is now constructed based on p. The set of agents in 
S' is 



the set of resources is /' = {i G / : Vi C V'}, and the set of beneficiary 
parties is K' = {k E K : V'}. The coefficients ajj^ and Ckv for 

i E I' , k & K' , and v E V are the same as in the instance §. The locally 
unique identifiers of the agents v E V are the same as in the instance 
S. (If we prefer globally unique identifiers which are a permutation of 
1, 2, . . . , |V^'|, we can add redundant variables to V.) 

4.4. The structure of Next we show that the structure of S' is 
tree-like, that is, there are no cycles in the hypergraph "K' defined by 
the instance S'; by construction, Oi.' is a subgraph of "K. 

For each q E Q, the subgraph induced by Tg in J{ is a hypertree. 
Furthermore, the subsets Tg form a partition of V. Therefore any cycle 
in "K and, therefore, any cycle in "K' must involve hyperedges which 



(3) 




v 
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cross between the subsets Tq and, finally, return back to the same 
subset. 

The only hyperedges which connect nodes in Tq and for distinct 
q,w E Q are the hyperedges of type III. There is at most one such 
hyperedge for any fixed q ^ w; this hyperedge corresponds to the 
edge {q, w} in the graph Q. Therefore a cycle in IK' implies a cycle 
in BQ{p,2r); this implies a cycle of length at most 4r + 1 in Q; by 
construction, no such cycle exists. 

4.5. A feasible solution of Next we show that there is a feasible 
solution X of S' with u = 1. Let u be the root node in 7p. By con- 
struction, u E Tp C V'. For each v G V, let x^, = 1 if d^'{u, v) is even; 
otherwise, let Xy = 0. See Figure [T] for an illustration. 

Because S' is tree-like, there is a unique path connecting u to v 
in "K' for each v E V. In particular, this path is a shortest path 
and has length d^/{u,v). It follows that each hyperedge in "K' has 
a unique node (that is, the node having the minimum distance to u) 
of its distance parity to u. Observe that hyperedges of resources and 
beneficiary parties alternate in paths from u. By the structure of S 
and §', it follows that the hyperedges of resources (type I) have a 
unique node with even distance to u. Therefore, ^„gy/ aiyX^ = 1 for 
each i G the solution is feasible. Analogously, the hyperedges of 
beneficiary parties (types II and III) have a unique node with odd 
distance to u. Therefore, ^^^y, CkvXv = 1 for each k G K\ implying 
u = l. 

4.6. The solution achieved by A in §'. Now we apply A to §'. The 
local radius-r view of the nodes f G Tp is identical in both § and S'. 
In particular, the deterministic local algorithm A must make the same 
choices x^ for f G Tp in both instances. 

As there is a feasible solution with = 1, the approximation al- 
gorithm A must choose a solution x with Yliv<^v' ^^v^v > l/a for 
all k G K'. 

We proceed in levels £ = 0,1,..., 2R — 1 of Tp. We study the total 
value assigned to the variables at level £, defined by 

s{^) = J2 ^- 

Recall that \Tp{£)\ = {dOy/^ for i even, and \Tp{£)\ = {dD^^-^y^d for 
£odd. 

Let us start with level i = 2R — 1, that is, the leaf nodes in 7p. 
For each v G Lp, there is a. k E K' such that = {v,f{v)} and 
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Ckv = Ckf{v) = 1- Therefore, by (|3j) and the fact that 5{p) > 0, 

1 1 d^D^~^ 

(4) S{2R -l)=J2xv = + - ^ (x, + > 

Next, we study the remaining odd levels £ = 2j — 1 for j = 1, 2, . . . , 
R-1. Consider the set Fp{2j - 1) = Tp{2j - 1) U Tp{2j). Observe that 
the hyperedges of type II occurring in Fp{2j — 1) form a partition of 
Fp{2j - 1). Each of the dW^'^ hypered ges in the partition has exactly 
one node in Tp(2j — 1) and exactly D nodes in Tp{2j). The coefficients 
Ckv of each beneficiary party k & K' associated with these hyperedges 
are 1/D for all v e V^. Thus, by the approximation ratio, we obtain 
the bound 

(5) S{2j - 1) + S{2j) = Yl c-^-^- ^ dW^/a. 

keK':VlCFp{2j-l) v^Vl 

Let us finally study the even levels I = 2j for j = 0, 1, 2, . . . , i? — 1. 
Observe that the hyperedges of type I occurring in Fp{2j) = Tp{2j) U 
Tp(2j + 1) partition Fp{2j). Each of the d^D^ hyperedges in the parti- 
tion has exactly one node in Tp[2j) and exactly d nodes in Tp{2j + 1). 
The coefficients of the resources i ^ I' associated with these hyper- 
edges are 1 for all v G V^. Thus, by the feasibility of x, we obtain the 
bound 

(6) 5(2j) + 5(2j + 1) = Yl 5^ a,,a:, < 

ieI':V{CFp(2j) vev{ 

Put together, we have, for j = 1, 2, . . . , i? — 1, 

(7) ^(1) < 5(0) + 5(l) < 1, 

(8) S{2j -1)1 dW^/a - S{2j) f S{2j + 1) - f 1 - dW^ 



a 



which, together with the assumption dD > 1, implies 

1 i S{1) i S{2R -1)- (i--\Y 

V j=i 



^RjjR-i _ 1 \ d^D^ - dD 



2a \ a J dD 



Therefore a>d/2 + l- 1/{2D) + {d + 2~ 2dD - l/D)/{2d^D^ - 2). 
Should we have a < d/2 + 1 — 1/{2D), we would obtain a contradiction 
by choosing a large enough R. This concludes the proof of Theorem [Tj 
The same proof with D = 1 gives the following corollary which shows 
inapproximability even if both a^^ G {0, 1} and £ {0, 1}- 
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Corollary 2. Let Aj > 2 be given. There is no local approxima- 
tion algorithm for Q with the approximation ratio less than Aj /2. 
This holds even if we make the following restrictions: G {0, 1}, 
Ck, e {0, 1}, = 2, Ai. = 1 and Af = 1. 

5. Approximability 

We have seen that the approximation ratio provided by the safe 
algorithm is within factor 2 of the best possible in general graphs; 
there is no local approximation scheme if Aj > 2 or A^ > 2. 

However, the graph in our construction is very particular: it is tree- 
like, and the number of nodes in a radius-r neighbourhood grows expo- 
nentially as the radius r increases. Such properties are hardly realistic 
in practical applications such as sensor networks; if nodes are embedded 
in a low-dimensional physical space, the length of each communication 
link is bounded by the limited range of the radio, and the distribution 
of the nodes and the network topology are not particularly patholog- 
ical, we expect that the number of nodes grows only polynomially as 
the radius r increases. We shall see that better approximation ratios 
may be achieved in such cases. 

Formally, we define the relative growth of neighbourhoods by 

\B:Kiv,r + l)\ 
7(r) = max — , ; — . 

We prove the following theorem. 

Theorem 3. For any R, there is a local approximation algorithm for 
([T]) with the approximation ratio 7(i? — 1) 7(-R) and local horizon Q{R). 

To illustrate this result, consider the case where !K is a (i-dimensional 
grid. In such a graph, we have \Bji{v,r) \ = 6(r'^) and \Bji{v,r -|- 1)| = 
|-Bj<:(f,T)| + B(r'^~^). Therefore 7(r) = 1 + 0(l/r) and our algorithm 
is a local approximation scheme in this family of graphs. 

We emphasise that the algorithm does not need to know any bound 
for 7(r). We can use the same algorithm in any graph. The algorithm 
achieves a good approximation ratio if such bounds happen to exist, 
and it still produces a feasible solution if such bounds do not exist. 
Furthermore, due to the local nature of the algorithm, if the graph 
fails to meet such bounds in a particular area, this only affects the 
optimality of the beneficiary parties that are close to this area. 

5.1. Algorithm. The algorithm is based on the idea of averaging local 
solutions of local LPs; similar ideas have been used in earlier work to 
derive distributed and local approximation algorithms for LPs (5l|7l|9]. 
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Figure 2. Definitions used in the algorithm. 

Fix a radius R = 1,2,...; the local horizon of the algorithm will be 
Q{R). For each agent u & V, define 

= B:Kiu, R), K'' = {keK:VkC \/"}, 

v;" = n 1/", r = {tei : ^ 0}. 

For each k & K and i ^ I, define 

Sk = f] V\ mu = \Skl Mk = max{|l^^| : j G \4}, 

= [jV\ N, = |[/,|, = min{|y^| : j G V,}. 

See Figure |2] for an illustration. For each u & V, let be an optimal 
solution of the following problem: 

maximise cu^ = min > CkvX^, 

k&K^ ^ 
v&Vk 

(9) subject to ^ a^^,x" < 1 for each i E 

> for each v G l^''. 

The solution x" can be computed by the agent u; or it can be computed 
separately by each agent j e which needs x", by using the same 
deterministic algorithm. 
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The agent j & V makes the following choice, which depends only on 
its radius 2R + 1 neighbourhood: 

(10) ft=mm|. = 

5.2. Constraints. Consider a resource i & I. We note that 

i G Vi and u <^=^ u eUi and j G Vi and u E 

(11) n G f/i and j G Vi and j G 

<^=^ M G t/j and j G 

and 

M G f/i ^ 3j G : ?i G ^ 3j eVf.jeV 

By definition, /3j < Ui/Ni for all i & Ij, that is, for all j G Vi. Combining 
these observations, we obtain 

Therefore x is a feasible solution of ([T]). 

5.3. Benefit. Let x* be an optimal solution of Q, with oj = u*. Then 
X IS cL feasible solution of (|9]), with co"^ > to*. Therefore the optimal 
solution of (IqI) satisfies 



(13) Yl ^'^-^ ^ ^* 

for all /c G -ft''^. Let /5 = min^gy jSj = mini^irii/Ni. 
Consider a beneficiary party k & K. We note that 

(14) u G Sk and j G Vfc =^ j G ^/c and u G 
and 

M G Sfc =^ u G K^' for all j G \4 j G for all j G 

(15) _ ..0, , ITsl 



C V^" ^ A; G ir" i=4 5^ Cfc^x^ > u*. 

vdVk 
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Combining these observations, we obtain 




In summary, the solution x approximates ([T| within the approxima- 
tion ratio maxkfzK Mk/mk-Taaxuzj Ni/iii. To complete the proof of The- 
orem |3| observe that max^^K Mk/rrik < 'j{R — 1) and maxii^j Ni/rii < 
7(i?). 
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